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ABSTRACT 


The results of an experimental investigation of the effect 
of adverse pressure gradient on the vortex breakdown phenom- 
enon in a diverging tube are presented. Adverse pressure 
“pradient was found to belas significant in determining the 
breakdown position as’ were the previously known parameters, 
namely, the Reynolds and circulation numbers. It was found 
that an increase in any of the three parameters serves to 
move the breakdown upstream, toward the origin of the vortex. 
It was further found that beyond some range, an increase in 
adverse pressure gradient causes boundary-layer separation 
in the diverging tube. Separation restricts the swirling 
flow and thus limits the effective adverse pressure gradient 
acting on the vortex. Data are presented to illustrate 
these effects. An approximate momentum analysis tor pre 
dicting the vortex breakdown position was carried out. The 
results were found to be only marginally satisfactory. In 
IW eeror thats however, the method appears to hold much 


promise as a breakdown predictor. 
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πι EN TRODUCT LON 


Vortex breakdown is an abrupt change in the NE 
of the core of a swirling flow. It was first noted by Peck- 
ham and Atkinson in 1957 during their studies of air flow 
over delta wings at large angles of attack [1]. Since then 
vortex breakdowns have also been observed in straight or 
Bua σας evlindrical ubes. Cyclone separators, behind Tor- 
pdecxondnpngborhnozzelscand dirrusepbns. Vortex breakdown 
ος ος tack OGeurean any vortex if the critical, swirl and 
flow conditions are met. Hence the study of this phenomenon 
may have a practical application in any device containing a 
Swirling fluid. 

Three distinct types of vortex breakdown have been re- 
ported by Sarpkaya [2]. These are the spiral, axisymmetric 
and double-helix types. Both the spiral and axisymmetric 
types initially form at a stagnation point along the vortex 
core, or centerline. These breakdowns mark the transition 
between laminar and ο... flows. 

In the case of Spiral breakdown, the core appears to 
Mr rsowstdadenliy stagnate, then kink COo one side, and then 
form a spiral as the kink rotates. After a few revolutions 
the spiral becomes lost in the violently turbulent fluid 
flow. 

The axisymmetric breakdown also initially appears as 


a stagnation in the flow, followed by a roughly right-angled 





MIO ENE Tr ather than forming espiral, the rotating 
kink forms the forward boundary of an egg-shaped body, which 
ie 115 from the rear, exhibits significant internal backflow 
i Ob lb Sor Ch ana then Gischarges through the®dewnstream 
Side of the body. Depending on the magnitudes of both the 
circulation and flow rate (i.e., Reynolds number) in the 
vortex the axisymmetric body may be closed from as little 
as 50% of its projected length,” up to 95%. When the body 
is only partially closed, turbulence commences at the exit 
of the body, and fans out to the outer reaches of the still 
laminar vortex in a very rough conical shape. When the 
axisymmetric "bubble" is almost fully closed, the exit core 
appears to re-form for one or two body lengths, and then 

επ ποτ Kink againwane! form a spiral as described above or 
proceed immediately to full turbulence in a smooth conical 
shape. 

The double-helix type of breakdown does not include a 
ΠΕ τ; POLLA  5ΠΕΝΟΟΠΕΙ rather, a mwd diis- 
turbance causes a slight slowing and swelling of the core. 
The core then shears into a tape-like helix formation. Its 
appearance is that of two ribbons emanating from a slight 


expansion of the core streakline, then entwining about each 


l This type is sometimes referred to as a "near- 
axismmetric" type [2]. 
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other. In the case of the double-helix breakdown, which 
generally occur at low Reynolds numbers (Re xz 700 to 1500), 
ΕΟ". 1.5 mediae transition from laminar to turbulent 
flow. It may therefore be considered by some to be more of 

a flow instability, significantly different in nature from 
the more common and relatively abrupt spiral and axisymmetric 
vortex breakdown forms [2]. 

Shortly after Peckham and Atkinson reported their dis- 
Over luumerous researchers attempted to develop a suitable 
model to describe and predict the occurence of vortex break- 
down. The experimental contributions of the various workers 
have been well documented by Robertson [3], and will not be 
ο ο σι ο... on Estas Chegezperimental ο ΤΟ, several 
theoretical analyses were developed, but the divergence of 
approach and initial assumptions have created significantly 
different and conflicting theories to explain the vortex 
breakdown phenomenon. Recently Hall [4] presented an excel- 
lent summary of these studies and compared the various anal- 
yses, highlighting their major strengths and weaknesses. 
Hewerouped the various) proposals into three categories, in 
which thes basic tldeas weretas follows: 

l. The phenomenon is in some sense like the separation 

of a two-dimensional boundary layer (Gartshore 
1962, Hall 1967, and others). 


o. The phenomenon is a consequence of hydrodynamic 
instability (Ludwieg 1962, 1965, and others). 


3. The phenomenon depends in an essential way on 


the existance of a critical state (Squire 1960, 
Benjamin 1962, 1967, Bossel 1967, 1969). The 


T T: 





Proposals ΠΠ thls "catlepory, havevequivatent 

ο μα OSO, Meter beca state but a 

Conce o μπε ρε Long Ot Usas "gen ticance. 
Hewewver, αἴνει Hmetalled fanedysis ob thestheories, Hall con- 
cluded that none of the proposals was adequate to completely 
explain what a vortex breakdown was, and when or where it 
would occur. 

Interestingly enough, Sarpkaya [5] noted that while the 
axisymmetric vortex breakdown basically corresponded to the 
finite-transition critical-state theory of Benjamin [6], 
the double-helix and spiral types were more in accord with 
Ludwieg's instability theory [7]. Carrying this one step 
further, Mager [8,9] presented a Em enatical analysis based 
on the works of Hall and Sarpkaya, which essentially com- 
bined the finite-transition and hydrodynamic instability 
theories. In essence he showed that each theory was basi- 
cally correct as far as it went, but that each was only a 
portion of a larger theory. He proposed that a complete 
analysis should, not stop when the breakdown occurred, but 
rather carry through from the laminar vortex to the onset of 
turbulent mixing [9]. This proposal will be discussed further. 

From the foregoing, it should be obvious that much is 
still not known about the vortex breakdown phenomenon. Not 
only is there controversy over how or why the breakdown 
Seecurs, but also on τας degree of importance of the various 
parameters; e.g., velocity profile, magnitude of swirl, Rey- 
nolds number and the effect of axial pressure gradients. 


Sarpkaya [2] has shown that the first three of these are 


Me 





quate Significant, at least in the study of breakdowns in 
water. Several researchers including Halil [H], Jones [10], 
Lambourne and Bryer [11], Mager [8] and Sarpkaya [2], have 
concluded that an adverse (positive) axial pressure gradient 
in the flow has an effect on the occurrence of breakdowns. 
Lowson [12], as a result of his analysis and comemiumenisss 
reported that "the final position of the majority of the 
practical breakdowns is due primarily to the reaction of 
pressure gradient effects on an already established break- 
down." However, as yet, no one has definitely verified either 
theoretically or experimentally the role that the axial 
pressure gradient is presumed to play. 

The purpose of the present investigation then was two- 
fold. The first objective was to examine what, if any, 
effect adverse pressure gradient EET vortex breakdowns. 
This was done by using cylindrical test sections of varying 
degrees of divergence to provide different pressure gradients 
Dos nes vor Leser eve ihe seeondzeb,jeer) verwasztresbr \elly 
examine Mager's theoretical analysis of vortex breakdown, 
and to determine whether it offered any more promise as a 
vortex breakdown predictor than: its predecessors. Considerable 
data were taken on vortex breakdown position and type over 
a wide range of flow ate Heb nd this was used to both 
Praphically demenstrate the effect of the ‘adverse pressure 
gradient and provide data for a computer analysis of Mager's 
hypothesis as suitably adapted to the present experimental 


system. 
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II. DESCRIPTION OF THE EXPERIMENTAL APPARATUS 


Bhewoecpenmmcnr-leopPapabuudsccensisteuRort a LESt tank, 
ot eS cio centrifuga punp Overhead reservoir, 
inlet box, two dye injection lines, vane-angle adjustment 
mechanism, thermometer, vortex-breakdown position/swirl 
pacucEmcEsmEICEdCVIcCEELwO potcameLeps amdgeneocessupy pups 
and valves, connected as shown in Figure 1. 

À scale drawing of the test tank is presented in Figure 
A was constructed of 3/4" Plexiglas plates, held to- 
gether by screws and sealed with silicon sealant. The re- 
movable top was made in two parts for ease of handling during 
Bono Vd. ΠΡ, {ΠΟΤ ee t Toan, and Was wot ted To 
mac ος cocoa asker ernal iuxw/lnes ari 


the test tank were as follows: 


l. Sectioned baffle-plate - This was made. of 3/4" 
Piles ilas ους πο σα εστω 59 rows ot' holes 


stepp d iron Pot avethe koccon to. 5S? 16" Bon 
top. This provided a nearly uniform velocity 
distribucion in tche test tCank: rine mesh 
screens were installed on each side of the 
οι ς plate" sections Merce dice chezrurbirrenee 
AL NENE Le a OA 


ο Adjustable vane mechanism An additional 3/4" 
paco ase pubes WoomobuaeciCOm be tne Upstream 
wall of the test tank. A movable circular ring, 
9" in diameter, was set into this plate. 
Thirty-two fixed pins were attached to the ring 
AE Sy miei cua Narra”. 11,25 degrees 
apart. Each of these pins fitted into a slot 
ΤΠ mp Lage. οσσαπε,.. Tne vanes 
Were” pinned to the inner plate; vand pivoted 
together about their fixed pins as the movable 
ring was rotated by a hand-screw mechanism. 

The indicator on the hand-screw mechanism was 
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calibrated to read the vane angle directly to 
the nearest half degree, and was capable of 
changing vane angle from plus to minus sixty 
degrees, thus varying the amount of swirl 
imparted to the vortex. Blade dimensions and 
geometry are shown in Figure 3. 


3. Cone and bell-mouth - A six inch diameter cone 
was attached to the additional plate described 
above, inside the circle of adjustable blades. 
This cone and its matching bell-mouth formed 
the entrance to the test section. The center- 
line was bored to provide an entrance for the core 
dye stream. An outer dye stream (to show the swirl 
angle) was introduced through an adjustable probe 
located in the flow stream between the core and 
bell-mouth. To ensure that the bell-mouth (and 
the connecting test section) were centered on 
the cone, simple hand-screw adjustments were 
Drevaded siGmmromim nor! Zemual and vertical 
centering. 


4. Divergent test section and outlet tube - Four 


pairs of these were made in matched sets, machined 
to close tolerances to ensure that there would 

be no obstruction to smooth flow through the bell- 
mouth, test section and outlet tube. To test 

Lho Er eC epressuremer adienu,estULDeCS or Varying 
angles of divergence were manufactured and in- 
stalled. The inlet diameter of all test sections 
was identical to that of the bell-mouth outlet, 
1.755". Test section #1 (TS-1) had a length 

of 10 3/8", and an outlet diameter of 2 1/4", 
providing a (half) angle of divergence (y) of 
1.38°. Test section #2 was made by adding a five 
inch extension to TS-1, maintaining the same angle 
of divergence, thereby giving it an outlet di- 
ameter or mr o T ned a length of 15". and 

v= ?.seen and, TS-ithed ag@lenctheor 15 3/01 

and y = 3.25°. Each test section was installed 
won “an outlet tube! that matched the Test sec- 
tion diameter. The outlet tube was mated with 

a@ matching nozzlé at its’ far end, to constrict 

the flow into the 1 1/2" outlet pipe leading 

to the rotameter bank. 


Exterior fittings on the test tank were the inlet box, 
thermometer to measure water temperature, and a 26" over- 


flow column to maintain constant tank pressure and absorb 


15 





flow fluctuations. A small device to measure breakdown 
position and swirl angle rolled along a section of model 
railroad track permanently attached to the top of the test 
tank. A one inch deep Plexiglas sighting bar was suspen- 
ded from a wheeled, counter-weighted cart. The sighting 
bar also provided a suitable location to mount a fixed pro- 
tractor and rotatable straightedge to measure swirl angle. 

To keep the system as free of disturbances as possible, 
air vents were provided at several locations, and the entire 
test tank was mounted on foam rubber. 

Flow of water through the system can readily be traced 
by referring to Figure 1. City water entered the lower 
reservoir, was pumped to the overhead reservoir, flowed by 
gravity into the inlet box, and entered the test tank. It 
then passed through the baffle-plate and screens, slowly 
traveled the length of the test tank, passed through the 
adjustable vanes, and between the cone and bell-mouth where 
ΙΙ πεα 5 Laminem νοχνζεχ. fhe water then entered the 
diverging test section and exited the tank through the outlet 
tube. After passing through the outlet nozzle, the water 
was carried to either (or both) of the rotameters, and 
discharged to a drain. Rotameter ratings at 100% flow were 
10.0 and 2.70 gallons per mane csneetl vel, and provided 
the capability of achieving Reynolds numbers within a range 
of Re = 300 to 18,000. 

A photograph of the system is presented in Figure 4. 


Out oft sieht are the pump, lower reservoir, inlet box and 
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the overhead reservoir. The large cement column to the 
rear of the system supports the overhead reservoir, located 


approximately fifteen feet above the test tank. 


ως 





III. EXPERIMENTAL PROCEDURE 


The initial phase of Tie investigation consisted of 
determining the location of vortex breakdowns in tubes of 
Various angles of divergence for identical flow conditions. 
ΠΕ ο πο μεσω οσοι ος Controlled were erreularıons 
Reynolds number, and degree of divergence of the test sec- 
tion. Circulation (2), a measure of the vortex strength, 
Was taken as the product of the tangential component of the 
velocity of the fluid leaving the vanes and the circumfer- 
ence of a circle drawn through the vane tips. From conti- 
nhuity and geometry, this resulted in the following equation 


[2], 


I = en(R, - m cos 6) V 


Lola s sin y 


where $ is the vane angle (read directly from the indicator 
in degrees) and the other terms are calculated or defined 
as shown in Figure 3. Circulation was normalized by defining 


a circulation number as 
Q = l/U. D. 


which may be expressed as 


1° D. (R_- m cos $) sin W 
we x 
64 adh I 


* 
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A simple computer program was written to calculate the 
corresponding circulation number for each vane angle 
eee rice 

Reynolds number (Re) was used to normalize the flow 


rate through the test section, 


ο ὃν 


ν 


Re = 





where D. was the test section diameter at the start of 
divergence (1.755"), U, was uniform velocity, and v the 
kinematic viscosity of the working fluid, water. U. was 
determined from the flow rate measured with a calibrated 
rotameter. As the experiment was run using city water, v 
did not remain constant, but varied as a function of water 
temperature. This variation has been taken into considera- 
tion in the calculation of Reynolds number by using the 
appropriate value of kinematic viscosity for the temperature 
measured. 

The experimental effort was conducted in four parts 
(one for each test section), and each part was conducted in 
an identical manner. The test tank was filled very slowly 
from the overhead reservoir in order to minimize turbulence 
in the test tank. Air vents, located at strategic points, 
were opened to eliminate any air pockets that might other- 
wise have caused disturbances in the breakdown position. 


After the test tank was filled, it was allowed to remain 
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idle for approximately one hour to permit the remaining 
Siam aneumunoie Turbulence To dissipate.* When the fluid 
was presumed still, the water supply was resumed and water 
run at a low flow rate through the entire system. 

Data runs were made by selecting a suitable vane angle 
and holding it constant while varying the flow rate through 
the test section by carefully setting the rotameter valve 
or valves. Vane angle setting varied from every two or 
three degrees for lower circulation numbers to five degrees 
forgniicgher circulations. The total range of circulation 
πο ΟΠ ΤΕΕ; Sections Varied from ? = Ὁ ο ο ο. 08 
($ = 13° to 55°), although only a portion of this range 
could be used for each test section. 

For a fixed vane angle, the rotameter setting was in- 
creased in iomall increments, and conditions in the test 
section recorded. When a vortex breakdown was present, the 
following information was recorded: i 
Vane angle (% - to nearest 0.2%) 

Flow rate (% - to nearest 0.5%) 

Water temperature (°F - to nearest 0.2°) 

Swirl angle’? (° - to nearest 2°) 

Breakdown type, or types if in hysteresis region 


Breakdown position -(inches from start of test 
section divergence - to nearest 0.1") 


ONO SW YF 


2 It would have been preferable to allow the water to 
remain in the test tank from day to day, but this would have 
permitted dissolved air to settle throughout the system and 
Gisturb the flow. 


3 Swirl angle was taken as the angle of the tangent to 
the outer dye stream at the nearest crossing of the core 
prior to the breakdown. It is considered by some to be an 
indication of a critical flow condition in the vortex [2], 
but it was inot used in this investigation. 
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Breakdown diameter (inches - to nearest 1/16") 
Percent that breakdown was closed, if axisym- 
metric type (% - to nearest 5%) 

Ratio of length to diameter, if axisymmetric 
type (% - to nearest 25%) 

10. Comments = anything peculiar to the current 
data wernt: 


O ON 


Breakdown and swirl visualization was accomplished by 
injecting dye into the vortex as previously described. 
Commercial food coloring at full strength was used — 
dye, and although it was known to be slightly denser than 
water, the effect off the density difference was not 
Delrsceable” 

Ric Os ıme een sumninge portion or the experiment was 
marking the location of the vortex breakdown. As the break- 
down is a dynamic phenomenon, its position oscillated within 
a narrow range, even though circulation and flow rate re- 
mained constant. All data points taken exhibited a slight 
ο ος cC puas lo spo Lie oXlamealreclion, he average 
distance traveled during these fluctuations was approximately 
tł 0.2 in. As experience was gained, it was found that by 
following the slow motion of the breakdown with the position 
indicator for about three minutes and mentally noting the 
limits, there would be roughly a 95% probability that the 
breakdown would remain in a 0.4" range. (This was true for 
both large and small flow rates.) The midpoint O πι 5 
range was taken to be the (average) breakdown position. 

This procedure was further complicated in the hysteresis 
region when two (or more) breakdown types were alternately 


present for the identical flow and circulation conditions. 
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When this occurred, it was necessary not only to determine 
the average position for each, but aiso the relative propor- 
tion of time that each was present. This latter step was 
necessary to facilitate the preparation of graphs in which 
average breakdown position for a given flow and circulation 
setting would be shown. 

Finally, after two complete sets of data runs were made, 
photographs were taken of various breakdowns over a wide 
range of flow conditions. These were taken on 2 1/4" x 
2 1/4" negatives in both color and black and white, using 
ο τι ποτε πως ασ to the rear of, and under the 
test tank. Several of these photographs have been reproduced 
P Pr il Oun #5, and $9 through 25. 

VROrseeonpietiomsor data runs and photography, vhe test 
tank was opened and cleaned, and a different test section 
fitted and installed. As the test sections were partially 
completed in advance, this step was accomplished within a 


day or two. 


Do 





IV. ANALYTICAL CONSIDERATIONS 


An exploratory approximate analysis of the effect of 
adverse pressure gradient on the location of vortex break- 
Πο consisted of the consideration of the equations o£ 
conservation of mass and momentum, the representation of 
Behe approximate velocity profiles by suitable polynomials, 
and the calculation of the distance at which the axial 
Velocity gradient increases rapidly leading to a stagnation 
της. 

The equations of motion for axially symmetric, incom- 
απ οσους TI OW, Wien Ppradiencs ın ché axial direction much 
πα) ο τισ ους τ πα τοςιαὶ αισοοῦσἹοῃ, may be wricven 


as [8] 


(ur), + (wr), = O (1) 
y? 

πο, πα P; (2) 
(uvr) + (vrô) = lr (w, - DI, (3) 
(uwr), + (wor), 9 -GDP, + weten), (4) 


where Re, = Q R,/v, and all velocities are normalized by 
the maximum velocity Q, all pressures normalized by the 


maximum dynamic pressure ρα. and all lengths normalized 
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Dy sheceradius) of the tube at the start of the diverging 


test section R... 


The above equations can be reduced to the following 


integral equations [8] 


IN CC (5) 

Lj! *& 1359! * IW! - (8/2)^*(a'/a) - 2 
where 

ó 
I, = f (W-w)r ar, 
O 
T = an - w)r dr 
D E 2 
(7) 
T = pw - ira 
πα u 
Š 2 
Ign = J/ v r dr, 
O 
and 
. A l 
DN (νο - WW!) + s) a’ — (Ap)' sz uc eee 
(8) 
with 
δ. 6 ᾶ 
ορ] δες E f — dr (9) 


ΡΕ = 
where a = 6 and S = Br 


and the normalized circulation respectively. 
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9 are the normalized core area 





The foregoing equations may be integrated further only 
if one assumes a set of suitable velocity profiles. This 
procedure is quite similar to that commonly used in employing 
the integral-momentum equation of Karnani One may do so 


by setting 


"Θὰ, ΟΡΟΙ, σας V(f, + 6f.), 


(10) 
V = F: 
where 
n= r/6, f. - n*(6 - n + 3n*) 
(11) 
f, - n(2 - nó) f = m(1 - η) 


These velocity profiles satisfy the boundary conditions at 

r= 0 and r = ó, and reduce to the swirling-potential-flow 
velocity distribution outside the vortex core. The parameters 
ο and B are initially unknown and will be determined later. 
Finally, to.complete the formulation of the analysis, one 
needs to express the conservation of massflow within the 


diverging test section by writing 


IA (12) 


το. an example of this procedure, see Schlicting [13]. 


> 





where A is given by 
2 
AS Clete, Can Yes (13) 


ZA is the axial distance, and y is the half-angle of test 
section divergence in Equation (13). Evaluating the inte- 
grals given by Equation (7) by the use of the velocity pro- 
files, Equations (10) and (11), and inserting the various 


T; j values in Equations (5,6,8,9, and 12), produces the 


following four equations: 


2 à 2 2 
[w (Kig + aK.. + K. 50 )} - [°/4alja' *[aW (K, 5 + 2aK, 5) Ja' 


I 


Ops 2 
+ [%T (0.4 + 28b. 4) J8" + [2aW(K, y + Κα + Κ ρα ) 


+ 0.1 aW(1 - a)JW' = [0] (14) 


[-0.2W(1 - a)]a' * [0.2aW]a' + [A - 0.2a(1 - a) ]W' 


= [-WA'] (15) 


[(r/a)°(1 ou b,58^) Ja! t [2ow* Ja! 
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* [-(b4, * 28b,,) T^/a]8' * [2a^w - 2u]w' 


= [48(1 - a)/(a Re)] | (16) 


[W(K5g + oK,, )]a' + [aWK4.]o' + [aW (b 24 + ob, ) Jg' 


21 
+ Ta(K,, + 0K,,)] = [2/Re] f (17) 
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The coefficients Kg and by j are eontalnedein Appendix 
A. 

The above set of equations were written into a computer 
Peorram ana solved simultUanecougslygtor at, at, 6', and W', 


ae οσα. BB, ana W by writing, tor example, 


= t 
Asay a; x Eb 


wich*theMinitial conditions 


This process was continued until a singular point was reached, 
forcing the coefficient matrix of Equations (14) through (17) 
omc Ome his condition occursiwhen a stagnation point is 
approached along the axis of the test section, and corresponds 
to that of vortex breakdown position [8,9]. 

tit should be®noted prior to a discussion of the utiliza- 
tion of the computer program that Mager attempted to calcu- 
reiche εεττ .:ν"ι:ε:-ονν!Ε τίν chrouehzcher use 
of a set of velocity profiles other than those used herein. 
An initial effort was made to use his formulas for the pre- 
diction of the vortex breakdown cion. but it was discov- 
ereditat thare were not only several errors in his final 


2 


fermutation, but? the velocity profiles chosen by him did 


2 Two typographical errors were found: (1) on page 
ΠΠ ΠΕΡΤΕ Weeecoefricient C „ Should be 2Wa/5 vice 2WA/5. 


(2) On page 651, Appendix, function f41 should inelude the 
number 2697, vice 26977. 
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Hoe Lead COmol point ofgsingularity even after the errors 
were corrected. Thus the analysis presented above, although 
Similar to that pioneered by Mager, differs in many respects 
andecercainly from thesobjectave and resuttis#discussed 
by Mager in reference [8]. 

The computer program writtenifrom thefabove analysis 
has been included in a separate section following Appendix 
B. It was written in FORTRAN (IV), and ean be most efficiently 
used on a time-sharing system communications terminal. 
Utilization of the program proceeded in the following man- 
Ii επ angele oi divergence for a particular test 
section was entered into the program. Then the program was 
compiled and executed. The program was written so that it 
would pause and request the initial values of a, B, 6, Re, 
ena, and the print spacing format; the latter term included 
e procedure" could be watched iteration by iteration, 
or by using a large enough number, would print only values 
at the location of core stagnation. 22 and Re were entered 
as the parameters of the particular breakdown position in 
question. The initial values of a and B were unknowns. 
They were found for this system by matching predicted break- 
down locations with those obtained experimentally for the 
same parameters. The procedure was to select a point, insert 
ehe corresponding o, Re, and 2, fixa at a= 1.0, and 
increase B until the maximum value of Z was obtained. Then 
Usineschzcse, plus the previous 6, Re, and Q; & was increased 


or decreased until the predicted value of Z corresponded to 
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that found experimentally. These computed values of a and 

B were used to predict breakdown locations for various values 
Of ò, Re, and R. As will be discussed later, the foregoing 
ΠΟΡΟΥ “παν: asi Capable of aceurately predieting the 
effect of the increasing adverse pressure gradient and 
Reynolds number. This was all that was expected from the 
analysis since no precise calculation could be carried out 
without knowing the initial velocity profiles. The measure- 


ment of the latter is a prohibitively complex problem. 
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V. PRESENTATION AND DISCUSSION OF RESULTS 


A. EXPERIMENTAL RESULTS 

Photographs of the three primary types of vortex break- 
downs, axisymmetric, spiral and double-helix, are presented 
in Figures 5 through 7. These breakdowns could be produced 
in each of the four test sections, although the occurrence 
of the double-helix type was much less common in the test 
sections with higher angles of divergence, TS-3 and TS-4. 

In addition to the three major types, several sub-types 
of breakdowns and core disturbances were also noted. The 
first of these to be discussed appeared as an instantaneous 
stagnation of the core flow. This occurred as a low Re/low 2 
phenomenon, and was quite unsteady. It simply appeared as 
a temporary slowing of the core that was immediately washed 
out before a breakdown could occur. As Reynolds number was 
¡Nerea sed Lince on pelo constant, its next type that 
occurred was an "open bubble," shown in Figure 8. It formed 
from a stagnation point and exhibited the initial kink found 
HmoOLMeGiemspiral and axisymmetric types: Iit differed, 
however, from both of these primary types as 1G exhibited 
some backilow behind the kink, but did not possess the energy 
required to fully close its forward wall as in the axisymmetric 
case DUC ο πως τισ του οΙροιἹστιος. the core dye washed 


out of the breakdown after one-half revolution of the ro- 


tating kink. A further increase in Reynolds number produced 
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the spiral and near-axisymmetric types in the hysteresis 
region. These phenomena are demonstrated in Figures 9(a) 
and 9(b). Figure 9(a) is a spiral type and 9(b) ae near 
axisymmetric type. These two photographs were taken at the 
same flow and vane angle settings; the photographs were 
taken about one minute apart. It should be noted that in 
the lower range of Reynolds numbers the spiral form pre- 
dominated, while the axisymmetric form was the more common 
and persistent of the pair for higher flow rates. 

The hysteresis effect appeared to be caused by the 
streamlined shape of the vortex breakdown, and the action 
in the test section during hysteresis can be shown by the 
following cyclical pattern. Assuming that the breakdown is 
first a spiral type, it will move slowly forward (toward the 
start of tube divergence) and backwards in the axial direc- 
tion. When it moves forward beyond the limits of its average 
position (about + 0.2"), backflow will suddenly start and 
chegaxısymmetcrie Lype wi Ὅουσπ. Τπε now streamlined shape 
will permit the an to move forward as the fluid 
accelerates around it: will §steady uin e new average 
position approximately 1" 5ο εἴ forward of that for the 
corresponding spiral type. As the near-axisymmetric type 
slowly expands and contracts, ice wil] respectively move 
Slowly forward and backward in the axial direction. POr 
some unexplainable reason, it will occasionally move back 
too far (out of its average range), and the breakdown will 


OB SS nenas tocecurs ωρα will contanuesto 
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rotate at the same speed, but the breakdown position in 
the test section will now permit only the formation of the 
spiral type. Then after a few moments of axial fluctuation, 
the spiral will move forward, re-form as an axisymmetric 
type, and the cycle will repeat. A similar hystenesis action 
was found to occur at very high Reynolds numbers Lor the 
near-axisymmetric and closed axisymmetric types. 

Still further increases in Reynolds number (to 
Re z 7000) produced the tightly closed axisymmetric type, 
shown previously as Figure 5. It has been deseribed quite 
comprehensively in Ref. 2, and thus a description of this 
phenomenon will not be persued further encina 

The above breakdown type transition pattern from 
instantaneous stagnation to the almost completely closed 
axisymmetric type was found to occur in all teskt sections. 
For the less divergent test sections (TS-1 and TS-2) at 
large circulation numbers, a completely different breakdown- 
transition pattern occurred. At extremely low Reynolds 
numbers (re z 800 to 1200) a core instability analogous to 
τ... stagnation occurred, but it did not INVOLVE 
a stagnation point. As shown in Figure ne core Etreak- 
line was deflected from the test section centerline, hea led 
slightly, then re-formed and once again entered the test 
section centerline. This offset portion ANTAS Core did 
not rotate, but rather remained stationary jin the position 
shown. This disturbance was accompanied by a tare increase 


in swirl angle, shown in the photograph. Notice that 
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turbulence did not occur in the test section downstream of 
the disturbance. It had originally been assumed that this 
disturbance was simply a double-helix type for which the 

dye preferred to illuminate only one tail: however, even 

when great quantities of dye were admitted to the core ae 
evidence of an obscured double-helix type was found. A 
double-helix type could be formed from the core deflection 

by increasing the Reynolds number. When this was done, the 
slight swelling of the core occurred roughly at the Point 

of the initial deflection described above. As in the core 
deflection case, the double-helix did not rotate - it appeared 
completely stationary. Although fluid was moving slowly 
through the test section, there was no visible fluid mecum. 
A further increase in Reynolds number at constant circulation 
caused the double-helix type to re-form into a rather flat 
near-axisymmetric type (Figure 11). This breakdown also 
appeared motionless, except for very slow fluid flow within 
the "bubble" area behind the stagnation point, and a slight 
wavering of the two tails emanating from the breakdown. 
a the breakdown did not rotate. A further increase 
in Reynolds number either caused a transition directly to a 
spiral or near axisymmetric breakdown, or in some instances 
at rather high circulation numbers (2 = 1.65 to 2.65), the 
formation of a rotating two-tailed axisymmetric typer 

Figures 12(a) and 12(b) show this type. The two photographs 
were taken seconds apart, and are shown to demonstrate the 


rotation of the tails. In this instance, the tails rotated 
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at approximately 39 rpm. Further increases in Reynolds number 
caused a transition to the near-axisymmetric type as noted 
above. Only very rarely was there sufficient circulation 
and flow strength to produce a completely elosed axisymmetric 
bubble for high circulation breakdowns in test sections with 
low angles of divergence. This is probably due to the fact 
natr τος high circulation and Reynolds number combinations, 
an axisymmetric breakdown occurs extremely close to the bell- 
mouth EA (forward of the start of test-section divergence), 
πιο tche vortex isi not completely formed. 

προ ορ ο ο6ακάσσαι σοσῃ. «Ποια in Pigure 135 > occurred 
lu E cavwiheEcaspoulactdsen numbers b appeared to be a rather 
ο ο Οσο. but rather than Tillinsg from Its 
closed downstream end, it filled from an axial, circumferencial 
τον Ιπεσαζεα approximavely vwwo-thirds of the way back from 
the stagnation point. The nearly stagnant fluid inside the 
body completed the closed downstream portion of the breakdown. 
It emptied from the ΕΙ ο πο ο ασ) σης rather chan 
from the rearmost point of the body as would normally have 
been expected. ligne Xe) lato ca Cm oii oe phenomenon Is 
Cine rods, 

Figures 14 through 17 show curves of the average 
normalized breakdown position as a function of Reynolds 
¡RAID Os lanl Leunlavronss torzeach or che Tour Vest 
sections. Each curve represents the average of two data 
τις ο οἱ πο tOo account forethe preference £or one of the 


two (or three) breakdown types in the hysteresis region. 
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Mia accion random curves have been selectedy and” for these 
all data points from which the average position curves were 
drawn are shown. These points demonstrate the scatter of 
the data for the two runs, the range and amount of hysteresis 
tor Giiferemt breakdown types, and finally the? transition 
through the breakdown sub-types as Reynolds number was 
changed. 
Figure 18 graphically demonstrates the significant 
er Oct of adverse pressure gradient on vortex breakdown. 
In this figure the curves show average normalized breakdown 
position for various Reynolds numbers for each of the four 
test sections for a single circulation number, Q = 0.67. 
The significance of the curves for TS-2 and TS-3 is obvious - 
as adverse pressure gradient increases, the breakdown occurs 
closer to the origin of divergence for the same Reynolds and 
circulation numbers. However, the meaning of the locations 
of the curves for TS-1 and TS-4 are not so obvious. From 
the above, it would have been expected that the curves for 
15-1 and TS-2 would coincide as they have the same angle of 
a that the curve for TS-4 would be far to the 
left of that for TS-3 since the angle of divergence is so 
ΠΠ realicen These vexpectavions obviously did not occur. 
Inasmuch as TS-1 was five inches shorter than TS-2, 
the transition from the adverse pressure gradient caused 
by tube divergence to the nearly zero pressure gradient 
of’ the following ‘straight section had the effect of slightly 


reducing the effective angle of divergence. Hence the 
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adverse pressure gradient acting on the breakdowns in TS-1 
was sede less than that acting on those of TS-2. 

The location of the curves for TS-! relative to those 
of TS-3 was at first quite baffling. It was only after 
the outer dye injector was used to illuminate the flow 
νεται. “Ες τος section boundary layer that the actual 
fluid action in the tube was discovered. Figure 19 shows 
that boundary layer separation was occuring in TS-4. Rather 
than the boundary layer being confined to a very thin layer 
of stagnant fluid along the tube wall, there was a much 
ΕΕ, Π.Ι Ἠνετπεεπιη: Cube wall and tche voartex. 
In this thicker boundary layer, reversed flow occurred, and 
cnis was actually the boundary layer that the vortex felt. 
Moeg tnne char is visibilfezgin the upper portion of tne 
Ecc ui cure 20.5nd 2l 13 the actual vortex boundary. and 
its angle of divergence is much less than that of the test 
ΤΟ dE e ia A ισα la o roja aa el 
τ ΙΕ; υειςς, De Tar the more notable as Is 
shows the early stages of the formation of a Taylor-Goertler 
πόδες within that boundary layer. 

Figure 24 shows graphically the effect of the adverse 
pressure gradient on vortex breakdown. Here average norma— 
lized positions of breakdowns are plotted as a TOI TION os 
the test section divergence for test sections 2, 3, and 4, 
for various Reynolds and circulation number combinations. 


NOV CMA > E ES IO WANDOS tons iñor TS 1 have been 


36 





included, and show that the effective angle of divergence 
for that test section was variable, but approximately equal 


to 1.30°, slightly less than its actual divergence, 1.38°. 


B. ANALYTICAL RESULTS 

The results of the exploratory analysis of vortex break- 
dad? tios presented im Figure 257 Tin this figures 
several’ computer predicted breakdown locations are plotted 
with those obtained experimentally. Although there is sub- 
Stantial deviation between predicted and experimental results, 
particularly at higher Reynolds numbers, the experimental 
curves are fairly well represented by the computer results. 
The primary reason for the deviations noted is attributed 
rostBezdifficulty of exactly reproducing the initial velocity 
profiles through the use of the parameters a, B, and ó. In 
add ea ma themattealtcomplextty of introdueing the 
boundary layer development along the test section wall and 
the difficulty of assuming an appropriate eddy viscosity for 
the turbulent vortex core further reduce the accuracy of 
this analysis. Taken together, the approximations made 
herein preclude a more precise analysis based on the quasi- 
steady state momentum equations. Thus this analysis can 
only be considered an attempt to explore the general influence 
of the adverse pressure gradient, rather than a model to 
accurately predict vortex breakdown location. 

It is presumed that the precise measurement of the initial 


velocity profiles by more sophisticated techniques (possibly 


S 





through the use of laser-Doppler annemometry) would signifi- 
cantly improve this analysis, and thereby produce much more 


accurate results. 


38 





VI. CONCLUSIONS 


J. sd0verse pressure gradient has a significant effect 
on vortex breakdown position. An increase in adverse pres- 
surezssracıent has the same eilect on breakdown position as 
an increase in vortex circulation or mean flow rate - that 
of shifting the breakdown location upstream, Toward the 
origin of the vortex. 

2. An increase in the angle of divergence of the test 
section beyond some specified point does not result in an 
increase in the impressed adverse pressure gradient.  In- 
stead, the boundary layer separates and reversed flow occurs 
τ Πε tube wall. This thicker boundary layer, vice the 
test section wall, then defines the outer extremity of the 
swirling flow, and thereby limits the effective adverse 
pressure gradient acting on the vortex. 

3. An increase in adverse pressure gradient does not 
directly affect the type of vortex breakdown created. It 
does, however, cause the constant circulation curves to 
shift down and to the left (toward the origin) on a plot 
of normalized breakdown position versus Reynolds number. 
curves sot cCherserrceulaysenenumbers are thus forced into 
ΓΕ τετ Reynolds number region o instability. This im 
turn restricts the formation of the relatively less stable 
Decsakoovnsseparticulariyschose of the double-helix and 


re la Led PCS. 
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y. In addition to the three primary vortex breakdown 
types (spiral, axisymmetric and double-helix), there exist 
numerous sub-types of breakdowns and core disturbances. 

5. The model developed to predict the location of 
vortex breakdowns performs only marginally well. The 
assumptions On which it 1s based, particularly that of the 
guit velocity profiles of the vortex, require further 
refinement. The method does, however, show promise, and 
Preovidesea direction for further investigation of this 


matter. 
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APPENDIX A 


COPPER ECIENS USEDE NEHE ANALYTICA SOLUTION 


K10 


= 0.1, k. = 
= 191/630, ke = 
= 23/840, ke = 
= 31/2520, - 
= kp, by 47 
= 2k¢, b,,= 
EE Dow 
- 0.5 -k - K, * ky, 
=k, - ky, 
= 2k, - k, - 0.5, 

bio ft By1B»> 
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111/315, kz = 1/3 
πι. ke = 7/6 
Bo p. Κο = 1/60 
3/20 , Kk, 5= 1/12 
2k y > by" Κα 
Ak] > DE 
ESI 
πμ τος ο 


151 2 
Koo” boo + b 418 
K57 0 





FIGURES 


APPENDIX B 


dund 





ALOAAISIU 





aunsseud e 





Ge 
$492 91109 04 
19|1π0 Ὁ u01322S 1591 
Jue] 1591 
uun | 02 
MO [44940 
ALOA43SƏ4 
peƏu4əA0 


|9A9| 4920M 


ΜΟ| 449Λ0 
snzeueddy ¡eyuawildadx3 ¿o Buimeug Sui] “| ϑωαη611 


H2 





szuLod 30A1Ld 
QURA JO sn30| 

BULA aLqeAou 
s}uod 94Ρ 





e - =“ < < oo" 4e e m omo 


k— v/e 4 


SALOAJISIA DÁP 





C OTT :9|/895) suoisuaurg pue xuej 3S9]  'Z 9Anb tJ 


A s x E 
| 149 >| 


PT 
(„38 = 43dap) 


xoq 39|U! ude 





191991 1 |. | 
ae UN [09 MOLJADAO iL n96 


35 





Figure 3. Vane Geometry and 





Dimensions 
E 20 h 
R, = 4.6875" 
ν 
> m = 1.6875" 
Ry d = 1.0" (depth of vane) 
@ = 5.625° 
total M 2g = 21R ,/32 (distance 
between pivot points) 
Vean h = 2gcos 4- 2m8 
a sin 
A e Fos $ - m/R, 
V | 
total 32 dh - 
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Photograph of Experimental Apparatus 


Figure 4. 
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Figure 9(a). Example of a Spiral Breakdown in the Region 
of Hysteresis (TS-1, y = 1.38°, Re = 2100, 
Q = 1.04) 





Figure 9(b). Example of an Axisymmetric Breakdown in the 
Region of Hysteresis (TS-1, y = 1.38°, 
Re = 2100, Q = 1.01) 
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Figure 12(a). Example of a "Two-Tailed Bubble" 
(TS-2, y -» 1.389, Re - 1300, Q - 2.66) 





Figure 12(b). Example of a "Two-Tailed Bubble" After 90° 
Of Rotation (TS—2, Y 1.38°, Re = 1300, 
Q = 2.66) 
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Figure 15. Breakdown Locations for Test Section #2 
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Figure 16. Breakdown Locations for Test Section #3 
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Figure 17. Breakdown Locations for Test Section #4 
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Figure 18. Constant Circulation Number Curves of 
Breakdown Locations for Test Sections 
#1 Through #4 
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Figure 25. Actual Vortex Breakdown Position Versus 
Computer Predicted Position 
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COMPUTER PROGRAM FOR PREDICTING VORTEX BREAKDOWN LOCATION 
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